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Techniques developed for handing inverse-power-law potentials in atomic physics are applied to the 
tensor one-pion exchange potential to determine the regions in which it can be treated perturbatively. 
In 5*-, P- and D- waves the critical values of the relative momentum are less than or of the order 
of 400 MeV. The RG is then used to determine the power counting for short-range interaction in 
the presence of this potential. In the P-and Z>-waves, where there are no low-energy bound or 
virtual states, these interactions have half-integer RG eigenvalues and are substantially promoted 
' ' relative to naive expectations. These results are independent of whether the tensor force is attractive 

, or repulsive. In the ^5*1 channel the leading term is relevant, but it is demoted by half an order 

• compared to the counting for the effective-range expansion with only a short-range potential. The 
' tensor force can be treated perturbatively in those _F-waves and above that do not couple to P- or 

■ I. INTRODUCTION 

'NT 
> , 

— , Since Weinberg 1] first proposed that the ideas of chiral perturbation theory (ChPT) could be applied to nuclear 

• forces, there has been a continuing debate over which parts of interaction can be treated perturbatively and which 
' can, or indeed must, be treated nonperturbatively. This has led to two widely-used schemes for constructing effective 

field theories (EFT's) to represent these forces.^ 

One is based on Weinberg's original suggestion P] and has been widely applied by van Kolck and collaborators 
[E^IS' It will be referred to here as the WvK scheme. In it one first expands potential using perturbative "Weinberg" 
power counting (like that in ChPT for mesons or single nucleons |^). Then one constructs the scattering amplitude by 
iterating the lowest-order terms: the leading, energy- independent contact interaction and one-pion exchange (OPE). 
I . The other scheme, developed by Kaplan, Savage and Wise (KSW) , starts from a nontrivial fixed point of the 
Q renormalisation group, which corresponds to a two-body system with an infinite scattering length. In the expansion 
^ , around this point, all pion-exchange forces as well as momentum- or energy-dependent contact interactions are treated 
^ ' as perturbations. 

At very low momenta, pion-range physics is not resolved and nuclear forces can be described just in terms of 
contact interactions. In this regime the two schemes, WvK and KSW, are equivalent, since iterating the leading 
contact interaction leads to the same power counting QjlailSli they just reproduce the effective-range expansion 
[lol[T]| . However at higher momenta, they treat OPE differently. There, problems with the perturbative KSW scheme 
have been identified by Fleming, Mehen and Stewart |l2j . In particular they have shown that the expansion is only 
slowly convergent in the ^Sq channel and, worse, it seems not to converge at all in the ^Si channel. This has led 
Beane et al. [Tjl to propose a hybrid approach, using KSW in the former and WvK in the latter. 

The strength of the OPE potential is given by the square of the pseudovector coupling which, to leadin g ch iral 
order, is given by = gj^m"^ / {IGir F^) . Empirical determinations lead to values for of about 0.075 [23. If 

we factor the nucleon mass out of the Hamiltonian, we find that the OPE potential contains the scale 

9aMn 

which is constructed from the nucleon mass and the pion decay constant. Since both of these are high-energy scales 
in ChPT, one would naturally take this to be a high-energy quantity. If A^r were much larger than m^, OPE could 
the be treated as weak for scattering momenta of the order of and the KSW scheme would apply. This would 
be the case in a world with much smaller up- and down-quark masses, where one would be much closer to the chiral 
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^ For reviews, see Refs. I^l^. 
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limit and m^r would be smaller than /jr, not just smaller than 47r/^ — the typical combination that appears in chiral 
expansions for processes involving at most one nucleon 0. Unfortunately, in the real world A^r is only about twice 
TOtt, and so we do not have a good separation of scales. This is what underlies the difficulties in building a useful 
perturbative EFT with pion-exchange forces. 

The first question is whether to iterate OPE or not. However, Nogga, Timmcrmans and van Kolck [T^ have shown 
numerically that iterated OPE cannot be consistently renormalised if the contact interactions are assigned the orders 
they would have in naive dimensional analysis. A similar observation is also made in footnote 5 of Ref. jl5l |. but 
without examining in detail the consequences for power counting. Hence, if we choose to iterate OPE, we are forced 
to address a second question: what power counting should we use for the resulting contact interactions? 

In this work, I address the first question by constructing exact solutions to the Schrodinger equation and examining 
them for nonanalytic dependence on the strength of the OPE potential. In the chiral limit, this potential has a tensor 
for m propo rtional to \/r^. Solutions can be obtained using techniques that have been developed in atomic physics 
UJ,LJ,UJ- As in the case of the simple l/r^ potential studied by Gao , the solutions in each partial wave become 
nonanalytic in the strength above some critical value for the dimensionless product of the momentum and coupling 
strength. This implies that the potential must be treated nonperturbatively in this region. 

For the in uncoupled partial waves, the critical values have been determined by Gao Here I apply these 
results to the tensor potential in coupled waves. Since the strength of OPE is given in terms of l/A^r, these results 
can be converted into a critical value for the relative momentum in each channel, above which OPE must be treated 
nonperturbatively. In channels involving waves with I < 2, these critical momenta are ^ 400 MeV, implying that 
OPE needs to be iterated in them, in agreement with the observations in Ref. In contrast, the critical momenta 
in channels that involve only waves with I > 3 are well above 1 GeV and so OPE can treated perturbatively. 

To answer the second question, I use the Wilsonian renormalisation group (RG) [T9| to determine the scale depen- 
dence of the interactions between two nucleons. In this approach, one imposes a floating cut-off A, lying between the 
low-energy scales of interest and high scales of the underlying interest. Demanding that the scattering amplitude be 
independent of the cut-off then leads to an RG equation for the effective short-range potential describing the physics 
that is not resolved at the scale A pol l2lll2^ . 

This cut-off could be imposed on a plane- wave basis but, in the presence of a known lon g-ra nge potential, it turns 
out to be more convenient to work in the basis of distorted waves (DW's) of that potential |22,|22|. For a long-range 
potential that is singular at the origin, the scaling behaviour, and hence the power counting, is controlled by the 
power-law dependence of the DW's for small r. In the case of the 1/r^ tensor potential, I find that the short-range 
interactions have half-integer RG eigenvalues (anomalous dimensions) and so their scaling is quite different from that 
given by naive dimensional analysis. These eigenvalues are all positive if the scattering is weak, but are smaller than 
they would be in the absence of the long-range potential. Compared to the usual power counting P|, these interactions 
are "promoted" to lower orders in the expansion in small scales and hence are more important for low-energy scattering 
than one would naively expect. This agrees with the numerical observations in Ref. |l4l |. Similar conclusions about 
the need for a modified power counting are drawn in Rcfs. 0,113,13' which use a related DW approach with a radial 
cut-off only. The RG analysis here provides the new power counting. It also shows that in waves where the tensor 
force can be treated perturbatively, the usual power counting still applies. 



II. SOLUTIONS IN THE CHIRAL LIMIT 



The long-range OPE potential has a central piece 

V.cir) = - fL.^— ('^1 • '^2)(ti • T2), (2) 

and a tensor piece 

V.Ar) = X ^ (3 + 3m^r + mlr^) 5i2(ti • T2), (3) 

where S12 — 3((Ti • f)((T2 ■ v) — cri ■ <J2- 

In the spin-singlet channels only the central, Yukawa potential V-j^c contributes. This has a l/r singularity at the 
origin. Even when iterated, this is not sufficient to affect the power counting for the short-range interactions |2lj . 
However, in the case of the spin-triplet channels, we have to deal with the tensor piece of OPE. This behaves like 
at short distances and so it is not obvious that it can ever be treated perturbatively. If we do iterate this interaction, 
the resulting nonpcrturbative short-distance physics can alter the power counting. 

In fact a singularity is sufficiently short-ranged that waves with low momenta do not resolve the singularity 
and their scattering can still be treated perturbatively. As discussed below, the critical value of the momentum above 
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which nonperturbative behaviour sets in is proportional to the scale Ajr in Eq. and increases rapidly with the 
orbital angular momentum of the wave. This is because in higher partial waves the centrifugal barrier "protects" 
low-energy waves from probing the singularity. 

At long distances the potential falls off exponentially as a result of the finite pion mass. This ensures that there are 
no nonanalytic terms in the scattering amplitude and so, at very low energies, the effective-range expansion [Tol [Tll | 
can be applied to it. However it does not alter the singular behaviour at sort distances and so the effects of that can 
still be analysed using the simpler form of the potential in the chiral limit {rriTr 0).^ 

In the chiral limit the tensor interaction has the 1 / 7-'^ form: 

y^Ar) = j^JsMr,.r.). (4) 

This and other singular inverse-power-law potentials also arise for systems in atomic physics, where techniques for 
solving the corresponding Schrodinger equations have developed."^ 

In applications in both nuclear and atomic physics these singular potentials should not be regarded "fundamental" . 
Instead they provide the long-range parts of effective interactions. At small separations the finite sizes and structures 
of nucleons or atoms cannot be ignored, and other, short-range interactions are important. The solutions to the 
Schrodinger equation with one of these potentials can be used as a DW basis for analysing the scale-dependence of 
the associated short-range terms. Provided they are viewed in this way, as pieces of effective theories, even attractive 
singular potentials are meaningful (contrary to the comment at the end of Ref. ^3-) 



A. Uncoupled channels 

The tensor interaction couples spin-triplet partial waves with Z = j ± 1, such as "^Si and ^Di, but not those with 
I = j. I consider the latter first. For these, the tensor operator is just = 2 and the potential has a simple 
form. The solutions to the corresponding Schrodinger equation can be constructed as series expansions in Bessel 
functions using the method of Refs. Ill El 03 ■'^ Solutions for the pure I /r^ potential have previously been obtained 
by Gao [H, but I recap some of the main features of the method here before applying it to the coupled waves. 

The radial Schrodinger equation describing the relative motion of the two particles has the form 

The strength of the potential is B3 = — 6/(Af„A^) in the uncoupled isospin-singlet waves, such as and ■^6*4, 
and B3 = 2/(MjvA,r) in the isospin-triplet ones, such as ^Pi and '^F^. It is convenient to re-express this equation in 
dimensionless form by introducing the coordinate x = pr and the coupling k — pM^B^, where p — ^/M^E is the 
on-shell relative momentum. It is also convenient to put the radial equation into a form similar to Bessel's equation 
by defining 4'{x) = x^/'^'il){x). After some rearrangement, the resulting equation is 

^{x) = - (p{x). (6) 

X 

In this form we see that the behaviour of the solutions is controlled by the single combination of the energy and the 
coupling strength, k. 

This equation can be solved analytically with the aid of the methods in Refs. [H. Ill Il8l|. by expanding 4'{x) in 
terms of Bessel functions as 

00 

(j){x) = ^ anJn+u{x). (7) 
n— — oo 



1 



dr^ 



2 d + 
r dr r-^ 



dx2 



dx 



^ If Att had been much smaller than m,r then this would not be possible; the exponential fall-off would mean that waves did not probe 
the singularity until their momenta were at least of order m„ . As discussed below, the critical momenta obtained in the chiral limit are 
greater than in all channels except the -^Si-^Di. Even in that the case, the critical momentum is of the order of m^/2 and so the 
conclusion that perturbation theory breaks down for momenta of the order of should still be valid, although the precise value of the 
momentum at which this ha ppen s will be different. 

^ For a recent review, see Ref. l26l . A review of older approaches can be found in Ref. l27l . 

* Other closely related methods exis t for solving equations with inverso-power-law potentials. These are based on Laurent |2flt l3(l IsH or 
Bessel-product expansions l32ll33l . 
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The shift in the order by v is needed because the interaction on the right-hand side generates secular perturbations 
that must be resummed (see also: Ref. [s^, Sec. 11.1). Substituting the expansion {Tj) into the equation ||HJ) leads 
to an infinite set of linear equations for the coefficients a„: 



2(n + iy + l) 



(n + vf - + i) 
By introducing 6„ = a„/(n + v) and 

!i{n + v) = 2(n + v) [n + vf 
the linear equations can be put into the more symmetric form 

K&n^i - fi{n + v)bn + nhn+i = 



2(n + z/- 1) 



-CXI < n < oo. 



oo < n < oo. 



(8) 



(9) 



(10) 



Following Refs. [13, 113 (^^^ ^1^° -R-^f- ^29]) one can solve these equations to get a representation of the ratios of 
coefficients in terms of a continued fraction. For positive n, I define the ratios 



Rn — 



bn 



Kb„ 



These satisfy the recurrence relation 



Rn — 



1 



fi{n + - K?Rn+i ' 
In terms of these the 6„ can all be related to bo and expressed as 



n 



\m—l 



In a similar way, for negative n I define 



and these satisfy 



Rn 



b-n 



Kb 



fl{~n + v)~ K^Rn+l 



The corresponding coefficients can be written 



b-n=[X{Rm] k"6o. 



\m—l 



Using these results in the remaining equation of Eq. H1U|I . with n = 0, gives 



fi{v)-K^ (Ri+Ri) 



0. 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



This is a nonlinear eigenvalue equation which determines the shift v. Alternatively, if one is simply interested in the 
value of v, one can look for the zeros of the infinite-dimensional Hill determinant of the coefficients in Eq. (|10|l (see 
Ref. [s^l, Sec. 7.5). If the set of equations is truncated to a finite number, this can be done straightforwardly with 
the aid of Mathematica [s^. Taking |n| ^ 20 is sufficient to determine the zeros to six significant figures, at least for 
small angular momenta (/ ^ 5). 

The resulting eigenvalue equation, in either version, is an even function of k and hence the roots of the equation 
are the same for both repulsive and attractive potentials of the same strength. From the fact that fiiv) is an odd 
function, it follows that 



Rn{-v) = -Rn{i^)- 



(18) 
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Hence if is a solution of Eq. H17|l . so is v. Also, the roots are periodic under addition of any integer to v. From now 
on, I shall use v to denote the root whose real part lies between I and ^ + 5- This ensures that coefficient is large 
in the expansion of the corresponding solution, 0'^+)(a;), at least for small k. A second, independent solution, (\)^~'^[x\ 
is obtained by replacing v by ~v. From the symmetry of the i?„ it follows that the coefficients in the two solutions 
are related by 



aL-j = (-iraW, 



if we choose 



(19) 



(20) 



The coefficient ao can be fixed by requiring that the solutions have the standard asymptotic normalisation of a Bessel 
function for large r. 

Solutions to these equations can be obtained using perturbation theory if k is small, corresponding to weak coupling 
or, equivalently, low energies. This is just the Born exp ansion, and it leads to expressions for the solutions as power 
series in k. The advantage of the Bessel expansion [iajll3)ll3 is that it determines the radius of convergence of this 
series, as well as providing analytic forms for the solutions in nonperturbative cases. 

For K = the order of the Bessel function is = ^ + i , and the solutions of the radial Schrodinger equation are 
just spherical Bessel functions ji{pr). Consider next very small values of k. For these we can approximate the ratios 
by their leading-order expressions from Eqs. (|11I15|I : 



Rn 



and the eigenvalue equation (fT7|l becomes 



1 



R„ 



= 0. 



fi{l + ,y) fii-l + iy) 



l + \ 



/i(^)-«' 

Since the solution v lies very close to Z + i in this limit, we can write 

5v. 

In partial waves with Z > 1, we can approximate fi{n + v) by 

fi{v) c -4(? + i)'<5i^, 
fi{n + jy) ~ 4n(/ + n+ i) (; + 1)) 

Using these in the eigenvalue equation we get, at order k^, 

A 3 
ov — — 



(21) 



(22) 



(23) 



16 (z-i) (z+i) {i + i)i(i + iy 

In the corresponding solutions to the Schrodinger equation, only the coefficent is large. 
For s-waves we need to be more careful since /o(— 1 + i^) is small in the weak-coupling limit: 

/o(-l + -6v. 

Using this in the eigenvalue equation, we find that the leading shift 5v is of order k, and is given by^ 



(24) 



(25) 



(26) 
(27) 



In the corresponding wave function, a_i is of the same order in k as oq. 



^ Compare Ref. |2(1 Sec. 2.5.2, noting that I have made a different choice from the multiple roots of the equation for v. For s-wave 
scatt ering by a pure potential, this dependence of the order of the Bessel functions on k leads to nonanalytic terms in the phase 

shift Il8ll36ll37l . However these are not present for OPE, even in the chiral limit, since the tensor potential vanishes in s-waves. 
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For larger values of k we need to iterate the recurrence relations for the i?„ but, provided k is small enough, we 
can still expand the results as power series in k. Similarly the solutions to the eigenvalue equation and hence the 
wave functions themselves can be found perturbatively as Born expansions in powers of k Since k cx P-B3, these 
expansions describe systems with either weak coupling or low energy. 

As |k| increases, these series converge more slowly and eventually a perturbative treatment becomes impossible. A 
definite upper bound on the value of n for which this occurs can be found by following the behaviour of the eigenvalues 
of Eq. (|17|) . As |k| increases, pairs of eigenvalues approach integer values from above and below until, at some critical 
value, they form degenerate pairs. Then, for larger values of \k\, they move off into the complex plane ^18] . The 
presence of a square-root branch point where this happens means that eigenvalue u cannot be expanded in powers 
of K above this critical value, Kc. If the eigenvalue cannot be expanded in this way then neither can the solutions to 
the set of linear equations (|10|l . and hence /t^ also provides an upper limit on the convergence of any perturbative 
expansion of the solutions. 

The values of Hc for low partial waves are listed in Table I. For large orbital angular momentum (/ ^ 20) Kc 
grows roughly as 0.3 Z'^. This can be understood if nonperturbative behaviour sets in when waves can penetrate the 
centrifugal barrier to radii where the and centrifugal potentials are roughly equal. To estimate where this occurs 
we can set 

E.i|±il.% ,28) 

This leads to 

~ + (29) 

which is consistent with the observed growth with 1. 



I tic 






0.318058 


1 


2.51811 


2 


8.33342 


3 


19.6983 


4 


38.6026 


5 


67.0469 



TABLE L Critical values of of the dimensionless coupling k for which the eigenvalues v form degenerate pairs. These agree 
with the values of esc in Table I of p^. where esc = k^/4. 



This dependence of Kc on the orbital angular momentum demonstrates that it is short-distance physics that leads 
to the breakdown of perturbation theory: the wave function has to penetrate the centrifugal barrier before it can 
"see" the singular core. The radius at which this happens is, from Eq. H28|l . of the order of 

where 



/33 = AUBs (31) 

is the length scale associated with the strength of the potential. These estimates of the critical value of k and the 
associated radius fail for small angular momenta. In particular Kc is nonzero for 1^0. This is because the 
potential is short-ranged in the sense that, apart from a single logarithmic term, its low-energy scattering amplitude 
can be expanded in powers of the energy |18.. 36. .37] (unlike the Coulomb or potentials). Hence a minimum 
momentum is required before its singular nature can be resolved. 

In the realistic case of finite pion mass, the long-range tail of potential is replaced by an exponential fall-off. Since 
the breakdown of perturbation theory is a short-range effect, it should not be altered qualitatively by a nonzero m^, 
especially in high partial waves. Provided the radius tq is much smaller than the finite mass will not even 

change the critical value of k significantly. In lower waves where vq comparable to or less than l/m,^, momenta of the 
order of are needed before the singularity is resolved and nonperturbative behaviour sets in. Hence in channels 
where the momentum corresponding to Kc is less than mjr, such nonpertubative effects are still expected for momenta 
of order TOtt, but not for much lower momenta. 
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The forms of the wave functions at short-distances will be an important ingredient in the RG analysis below and so 
I outline their basic features here. More details of the solutions can be found in Ref. ^18\, for the case of repulsive 
potentials. For small r, we can use the WKB approximation to find their forms. In th e case of a repulsive potential 
this shows that the solutions have exponential dependence on -y/ k/x, or equivalently \f^zjr. In dimensionless form, 
the small-x solutions are 



exp 



X 



B^^\k)x^I^ exp 



-2,/- 

X 



(32) 



The determination of the coefficient A of the dominant piece of a solution requires careful asymptotic analysis of the 
series Eq. ((TJ. The important terms at small x are those for large negative values of n. These can be summed using 
Laplace's method (see Ref. [S^l, Sec. 6.7). Under the analytic continuation x e^'^^x the dominant and subdominant 
pieces of the solution exchange roles (an example of Stokes' phenomenon — see Ref. |34|, Sec. 3.7). This shows that 
their coefficients are related by 



(33) 



Having found the small- a; forms of the two independent solutions (x), we can then build the regular solution 
(which behaves like x^/'^ exp[— 21^^/2;]) as a linear combination of them. These regular solutions form a complete, 
orthogonal set of basis functions for the RG analysis of the short-range physics. 

In the case of an attractive inverse-cube potential, we need to be more careful. Since, as already noted, v and the 
ratios of determinants are independent of the sign of the potential, solutions for this case can be obtained by replacing 
K by — K in the recursion relations for the ratios of coefficients a„. Alternatively one can make an analytic continuation 
X of the solutions already found for the repulsive case. The latter method shows that the solutions have the 

small- a; forms 



cos 



2W- 

X 



1 

±1/- - 
4 



where I have now defined the dimensionless coupling to be positive-definite: 
shows that the C(k) are related by 

cW=2AW, 



(34) 



— pMjvl-Bal. The continuation also 



(35) 



to the coefficients A{k) of the dominant pieces of solutions with the same v for a repulsive potential of the same 
strength. 

Both of the solutions for an attractive potential display oscillatory behaviour as a; — *■ and so any linear combination 
of them is an equally good solution. For the RG analysis we need a well-defined set of orthogonal basis functions. 
As in the case of an attractive inverse-square potential,^ this can be obtained if we choose a self-adjoint extension of 
the original Hamiltonian. In practice this means fixing the phase of these short-distance oscillations ^2^,^3§\. This 
phase should be independent of energy to form an extension whose eigenfunctions are orthogonal. In essence one 
short-distance parameter, the leading, energy-independent term of the effective potential, has been used to provide 
a well-defined set of DW's of the long-range potential. Clearly this leads to a redundancy in the parametrisation: 
a different choice of extension can be compensated by changing the leading term in the short-range potential 2^. 
However any energy- or momentum-dependence due to short-range physics can be described entirely by higher-order 
terms in that potential. Note that here I am requiring orthogonality simply in order to generate a suitable DW basis for 
studying the effects of the short-range interactions. This is in contrast to the approach developed in Refs. 
where an orthogonality condition is imposed on the full wave functions and this leads to very strong constraints on 
the short-distance interactions. 

The forms of the short-distance wave functions for attractive and repulsive potentials look very different, depending 
sinusoidal or exponentially on y^fS^/r. However we shall see that it is their power-law radial dependence that controls 
RG flow of short-distance interactions and this is same for both cases. Although the wave functions for the attractive 
potential oscillate at short distances, these oscillations depend on a scale, P^, unlike the analogous ones found for the 
inverse-square potential. The scale-free oscillations found there and in the corresponding three-body systems lead to 
limit cycles in the RG flow 22, 38, 39, 40, 41, 42,^], but here the scale-dependence of the potential means that 
we should not expect to find similar limit cycles. 



See Refs. I22ll39l and references therein for more discussion of this potential. 
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B. Coupled channels 

The solution of the Schrodinger equation for the coupled spin-triplet channels proceeds along very similar lines. 
Using the matrix elements of in the two-component basis of waves with I = j ± 1 (2^, the chiral limit of the tensor 
potential can be written in the form 

""^^ ^ 2j + l[6VWTT) -20 + 2) ) r3' ^'"^ 

where Bt = — 3/(MjvA7r) for isospin-singlet waves, with I = j±l even, and Bt = 1/{MmXti) for isospin triplets, with 
I = j ± 1 odd. Rescaling the equation as above, it can be written in a similar dimensionless form to Eq. ©: 



' 2 d n' 

X —J + X \- X 

ax^ ax 



0-(L, + il)'0 = iK,0, (37) 



where the 2x2 matrices are 



' \ j + lj' ' + 6v/7(7TT) -2(j + 2) y ' ^ ^ 

and the dimensionless combination of momentum and coupling strength is 



-Sp/Ajr isospin singlet 
+p/XTr isospin triplet 



(39) 



The solutions to these equations can be expanded in Bessel functions as 

oo 

4>{x) = ^ anJ„+u{x), (40) 

n— — oo 

where the a„ are two-component vectors. Substituting this into Ea. (|37|l leads to an infinite set of linear equations. 
As in the uncoupled case, these can be put into a symmetric form by defining b„ = a„/(n -I- v) and 

Y{n + v)^['' ^ . (41) 

V f,+^{n + v)) 

The resulting equations are then 

Kb„_i - F(n + !^)b„ + Kb„+i = -oo<n<oo. (42) 

If we write 

b„ = R„Kb„_i, (43) 

then the satisfy the recurrence relation 

R„ = [F(n + KR„_iK]-i . (44) 
In terms of these matrices, the coefficients b„ for n> 1 are 

b„ = ( n R™K J bo, (45) 

\m=l / 

where the matrix product should be read as starting with m = 1 at the right. Similarly we can write the coefficients 
b_„ as 



W R„.K bo, (46) 
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where the R„ are given by 

r:„ = [F{-n + KR„_iK 
Consistency of the equation for n — requires that v satisfy 

det \F{iy) - K f Ri + Ri) K 



0. 



(47) 



(48) 



This is equivalent to the Hill determinant corresponding to Eq. H42|l . As in for the simple inverse-cube potential, the 
roots of this equation are symmetric in v and periodic. 

The 2x2 nature of the problem means that this eigenvalue equation has two roots vi and 1^2 with real parts between 
I and / + 5. The solutions to the Schrodinger equation (|S7jl for these roots will be denoted by cf)''^'~^\x) and cf)^'^'^\x), 
respectively. There are also two other independent solutions (f)''^'~\x) and 4>'--^'~\x), corresponding to the roots — I'l 
and —h'2- The starting coefficients bg ' must be eigenvectors of the n — equation 



F-K(^Ri + Ri)k 



= 0. 



Since F{n + v) is an odd function, the matrices for positive and negative roots are related by 

Rn(-l') = -R„(j^). 

If we choose the starting coefficients such that 



a, 



(i,-) _ Ji.+) ^ Ji) 



^0 — ^0 

then the coefficients in the pairs of solutions are related by 



= a, 



(-l)"a; 



' 



(49) 



(50) 



(51) 



(52) 



The aQ can again be fixed using the asymptotic normalisation of the solutions. 

As for the case of uncoupled channels, the eigenvalues v move into the complex plane for large enough values of 
and the solutions can no longer be expanded perturbatively. The critical values for which this happens are listed in 
Table II. The coupled nature of the equations means that there are in general two of these values. The one exception 
is j — 0, where only Z = 1 is possible and the critical value of is just a quarter of the one in Table I for a P-wave. 



j 


I 









1 


0.629528 




1 


0,2 


0.683495 


2.48290 


2 


1,3 


1.61857 


6.91983 


3 


2,4 


3.95647 


14.3624 


4 


3,5 


8.02206 


23.2001 



TABLE II: Critical values of of the dimensionless coupling kt for which eigenvalues form degenerate pairs. 

The same methods outlined above can be used to find the coefficient vectors a„ and hence to construct solutions 
to the radial Schrodinger equation (|37|l . Again, the full forms of these are not needed here, just their short-distance 
behaviours. These can be obtained by applying the WKB method to the eigenchannels of the potential. Since the 
eigenvalues of Kj are +2kt and —4:Kt, one of these channels is repulsive and one attractive. For small x, all the 
solutions are mixtures of exponential and sinusoidal pieces, with the form 



cj){x) ~ e+ ^(kt) a;^^'* exp 

+ e_ C{kt) x^I^ cos 
where and — k_ are the eigenvalues of K^, 



+ e+ B{kt) x^/"* exp 





Qp/X-n isospin singlet 
Ap/Xt^ isospin triplet 



(53) 



(54) 
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and e± are the corresponding eigenvectors. These eigenvectors are 

for an isospin triplet channel, and vice versa for a singlet. As in the single-channel case, analytic continuation in a; 
can be used to relate B(kt) to A{kt), and to express C{kt) in terms of the A{kt) of the solution for a repulsive 
potential. 

A basis set of physical, orthogonal solutions can formed by taking two linear combinations of these four independent 
solutions at each energy. These must be re gula r, with no admixture of the divergent x^^^ exp [2 / x] piece in the 
repulsive channel, as discussed in Refs. p^ l24| for the ^Si-^Di waves. As for the attractive inverse-cube potential, 
one energy-independent parameter is also needed to fix the phase of the oscillations in the attractive channel |l3j| . 
Although these solutions contain both oscillatory and exponentially decreasing pieces at small x, what matters in the 
RG analysis will be their power-law behaviour, which is the same as for the uncoupled l/r^ case above. 



C. Critical momenta 



For the one-pion-exchange potential of interest here, the strengths of the tensor interaction in the various channels 
are all fixed in terms of the scale A^r introduced in Eq. (Q. The critical values of the dimensionless coupling can 
then be converted into critical values of the relative momentum for each scattering channel. Taking — 290 MeV 
leads to the critical momenta listed in Table III. These are the maximum values of the momenta for which one could 
attempt to construct a perturbative expansion of the solutions. In practice one would expect such an expansion to 
be sufficiently convergent to be useful only for momenta well below these values. 



Channel 




Pc 




66 


MeV 


'Po 


182 


MeV 


'Pi 


365 


MeV 


''P2-F2 


470 


MeV 


'D2 


403 


MeV 




382 


MeV 


'Fs 


2860 


MeV 


'Fi-'^Hi 


2330 


MeV 




1870 


MeV 



TABLE III: Critical values of the relative momentum at which pairs of eigenvalues become degenerate and hence the tensor 
potential cannot be treated perturbatively. 

From the table we see that in two channels, ^Si-^Di and ^Pq, the critical values are less than or of the order of m^r. 
These are obtained from the chiral limit of the OPE; the values for finite ttItt will be somewhat higher. However they 
will still be of the order of m-^ since, as discussed above, the exponential fall-off cannot make the radius at which the 
waves probe the singular core much smaller than l/m^. It is thus not surprising that Fleming, Mehen and Stewart 
[T^ found that the perturbative KSW approach fails for these cases. In the other P- and I?-waves the values are low 
enough, ~ 400 MeV, to suggest that these should also be treated nonperturbatively at the energies of interest for 
nuclear physics. Again this is in accord with the findings of ref. [9]. 

Between the D- and P-waves, the critical momenta jump by a factor of about 6. Two effects contribute to this: 
the critical dimensionless couplings are about twice as large for the F-waves, and their physical couplings are three 
times smaller since the F-waves are isospin triplets. Hence for the F-waves and above that do not couple to P- or 
D-waves, the breakdown scales are well above 1 GeV. In these cases it should be possible to treat OPE perturbatively. 
A similar conclusion is reached in Ref. , but on the basis of very different arguments which rely heavily on keeping 
TOjr finite. 



III. RENORMALISATION-GROUP ANALYSIS 



The long-range, pion-exchange physics in low-energy nuclear EFT's can be calculated from ChPT 

Qai- The 

short-range physics is then parametrised in terms of contact interactions. If the EFT is to have any predictive power, 
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we must be able to organise these interactions systematically according to some power counting. In weakly interacting 
systems one can do this by simply counting powers of the low-energy scales, generically denoted here by Q. This 
naive dimensional analysis leads to the power counting used in ChPT for mesons and single nucleons 0. In contrast, 
nonperturbative effects can introduce new low-energy scales and can generate anomalous dimensions for terms in the 
short-distance potential. As a result these terms may not scale as predicted by naive dimensional analysis and so the 
power counting can be quite different. 

The RG provides a general and powerful tool for analysing scale dependences, particularly in nonperturbative 
systems. In the context of few-body systems, it is convenient to express this in the form of a differential equation 
describing the "flow" of the short-range potential as the cut-off is varied [23,0, 112 • This equation can be constructed 
as follows. First, we apply a floating cut-off at some scale A that lies between the low-energy scales of interest and the 
scale of the underlying physics Aq. This assumes that these scales are well separated; if they are not, the expansion 
in powers of Q/Aq will not converge and we shall not be able to construct a useful EFT. Second, we demand that 
physical observables be independent of the cut-off A, since its value is arbitrary. As a result, the couplings in our 
effective potential must depend on A. This dependence ensures, in particular, that the couplings cancel any parts of 
loop integrals that diverge for large A. Finally, we rescale the theory by expressing all dimensioned quantities in units 
of A. Powers of A can then be used to determine the net powers of low-energy scales in the terms in the potential. 

Having constructed the RG equation for the rescaled potential, we can look for fixed points. These are A-independent 
solutions that can form the end points of the RG flow. They describe scale-free systems. Perturbations around a flxed 
point can be expanded in eigenfunctions of the linearised RG equation that scale with definite powers of A, given 
by the eigenvalues of this equation. These perturbations can be classified according to their eigenvalues as relevant, 
irrelevant or marginal. 

Marginal terms ( "renormalisable" ones in field-theoretic language) have no power-law dependence on the cut-off 
after rescaling, although in general they can depend logarithmically on it. Like the fixed-point potential, these terms 
are important at all scales. They are the leading-order terms in the WvK scheme, of order 1/Q in low-energy scales. 
Since the loop integrals in the Lippmann-Schwinger equation are of order Q, all iterations of these terms are of the 
same order and hence they need to be treated nonperturbatively. Irrelevant (or "nonrenormalisable" ) terms vanish 
as positive powers of A as A —> 0. These are higher-order terms which become weak at low energies and so can be 
treated perturbatively. Lastly, relevant (or "super-renormalisable" ) terms grow as negative powers of A. These are 
unimportant at high energies but become increasingly important at low energies, ultimately changing the nature of 
the low-energy EFT. If such terms are present, a fixed point is unstable and, for low-enough values of A, the theory 
will ultimately fiow to a different point. An example of this is provided by the scattering- length term in the pionless 
EFT 2^. 

The RG approach developed in Ref. pl| assumes that two-body potential consists of a known long-range piece 
Vl and a short-piece Vs which parametrises the physics that lies outside the scope of our effective theory. This 
method starts by using the "two-potential trick" [ij to define a T-matrix describing scattering between distorted 
waves (DW's) of the long-range potential. The Hilbert space is then reduced by imposing a cut-off on the basis of 
DW's at momentum A. For this cut-off to lead a well-defined space, the DW's should form a complete, orthogonal 
set of basis functions. This requires that the long-range potential, and any parameter needed to form a self-adjoint 
extension of it, should be independent of energy. 

Demanding that the (fully off-shell) T-matrix be independent of the cut-off A leads to the differential equation 

where Gl is the DW Green's function for the long-range potential. In order to treat cases where the DW's vanish or 
diverge at the origin as a result of nonperturbative effects of , the short-range potential is taken to have the (5-shell 
form 

ys(p,A,A,i?;r) = l-5(p,A,A,i?)^^-^, (57) 

where A denotes a generic low-energy scale associated with V^. The radius R provides a second regulator here, which 
can be thought of as a "factorisation scale" separating the long-range physics which lies within the domain of our EFT 
from the unknown nonperturbative physics at shorter distances. This second scale is introduced because it leads to an 
RG equation from which the scaling of the potential can be deduced in a particularly transparent way. In approaches 
based on a simple momentum cut-off, as in Ref. or a coordinate-space regulator, as in Refs. 15, 23, 24], a single 
scale plays both roles of renormalisation and factorisation.^ For present purposes, R should be chosen to be small 



This is analogous to the choice that is often made in studies of the QCD evolution of structure functions. 
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enough that the wave functions have reached a common, energy-independent form. 
With this (5-shell form for Vs, Eq. (|56|) becomes 

^ - - ^ ^siP, A, A, R), (58) 

where ipLip,r) are the DW's for Vl and I have assumed that it does not produce any bound states. Taking the 
"factorisation" radius R to be small enough that it lies in the asymptotic region, the DW's can all be written in the 
form 

V-Lb, R) - N{\/p){pR)'^^-^^'^F{\R). (59) 

This is a slight generalisation of the cases considered in Ref. [23| to potentials that generate a non-power-law depen- 
dence on R in this region, described by the function F{\R). 

Following the general method of Ref. |2lj |. I introduce a rescaled on-shell momentum, p — p/ K. Other low-energy 
variables are treated similarly, and I then define the rescaled potential, 

Vs(p, A, A) - ^ (Ai?)--i \F{K\R)\^Vs{^p,^\K,R)■ (60) 

Inserting this and Eq. H59|l into the differential equation Eq. H58|l gives the RG equation for V^: 

A -r-p = p ^r:r + A — — + (tVs + — • (61) 

dK dp dX ^-P^ 

Since the function |F(AAi?)p depends on the product AA, its derivatives cancel in this equation. The evolution of Vs 
thus depends only on the power-law part of the DW's. 

Note that in Eq. (|60|l I have implicitly demanded Vs depend on R in such a way that the rescaled potential be 
independent of R. This implies that |V'l(p, -R)pVs is also independent of i?, and hence it ensures that scattering 
observables do not depend on this arbitrary radius. This is very similar to the philosophy adopted by Pavon Valder- 
rama and Ruiz Arriola in Refs. 0,|23,l3|i except that there the coordinate-space regulator is the only one. However 
the complicated dependence of F{\R) on R means that the scaling behaviour of the terms in the potential cannot 
easily be determined from this condition. In the present approach the radial regulator separates off the regime of 
nonperturbative short-distance physics, while the momentum cut-off A is used to analyse the scale dependence. For 
this to work, R must be in the region where the DW's iPl{p,R) have reached their energy-independent asymptotic 
form, at least for p < A. 

The starting points for analysing the scale dependence of the potential are the fixed points of the RG, solutions to 
eq. (I61() that are independent of A. Expanding Vs around one of these points and keeping only linear terms in RG 
equation, we get an equation whose eigenfunctions scale with definite powers of A. Since the rescaling means that the 
power of A counts the net power of all low-energy scales in each term, these eigenvalues determine the power counting 
for the terms in the potential. For example, in the case of a pure short-range potential 23|, the power counting 
for the expansion around the trivial fixed point is just the one originally proposed by WeinbergP|. There is also a 
nontrivial fixed point, which describes a system with an infinite scattering length. The expansion around this point 
can be organised according to the power counting developed in Refs. The terms in this are in one-to-one 

correspondence with the terms in the effective-range expansion [lol ITU . 

The form of the RG equation above can easily be used to analyse scale dependences in the vicinity of the trivial 
fixed point, Vs = 0. To find and study nontrivial fixed points, it is more convenient to convert it into a linear equation 
for l/Vs Mil- 

d\\vj apKvsJ a\\vs/ vs i-p' 

The detailed forms of the fixed-point solutions to this equation are not necessary to determine the power countings 
for perturbations around them. If needed, they can be found by applying the methods of Refs. |2H23)I3- These rely 
on the fact that the basic loop integral from the Lippmann-Schwinger equation satisfies the A-independent version of 
Eq. (|62|l . This integral contains a piece that is a nonanalytic function of p/X. Cancelling off this piece then leaves a 
well-behaved solution, Vso, which is analytic in the low-energy variables, p and A. 

So far I have discussed only energy-dependent perturbations. More generally, as an off-shell quantity, the short- 
range effective potential can also depend on momenta. This is equivalent to including terms with spatial derivatives 
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in Vs- For the expansion of a pure short-range potential around the trivial fixed point, the power counting is the 
same for both energy- and momentum-dependent terms. Hence one can "use the equation of motion" , to exchange 
one dependence for the other |45j |. 

The expansion around around a nontrivial fixed point is more complicated. There the equation of motion involves 
the fixed-point potential (which can include both long- and short-range pieces) and so purely momentum-dependent 
terms are not eigenfunctions of the linearised RG equation, in contrast to the energy-dependent ones discussed above. 
The coefficients of the latter terms are directly related to on-shell scattering observables through an effective-range 
expansion, as discussed in Refs. ^2(1, .21j . There are also terms with mixed momentum- and energy-dependence 
that scale with definite powers of A. However these affect only the off-shell behaviour of the scattering amplitude. 
Furthermore, they are of higher order than the corresponding energy-dependent terms |20|| . Hence, if one expands 
a momentum-dependent term in eigenfunctions of the RG, the scaling with A of the dominant piece is governed 
by the eigenvalue of the corresponding energy-dependent term. This is w hy t he analyses of Refs. QiHI^li using 
momentum-dependent potentials, arrive at the same power counting as Ref. [20|, which uses energy-dependent ones. 



A. Central OPE 



The first step in constructing any EFT is to identify all the important low-momentum scales. To illustrate the 
choices involved in the nucleon-nucleon system, I consider first the central piece of the OPE potential and summarise 
the relevant results from Ref. For scattering at energies ~ 100 MeV the relative momentum and the pion mass 
obviously form two of these scales. For the central Yukawa potential, we can construct the scale = m^/Ajr — 
70 MeV. In strict chiral power counting this would be of order Q^, since it contains two powers of m^r. However if we 
choose to treat as an additional low-energy scale, then is promoted to order Q. 

To see the consequences of this in the context of the RG for the ^5*0 channel, we can multiply the effective 
Hamiltonian by Mjv and define the dimensionless potential 

Then we need to express all low-energy scales in units of A. If we regard A^r as a high-energy scale, then we define a 
rescaled on-shell momentum hy p — p/A, and similarly TO^r = TOtt/A and f = Ar.^ The resulting rescaled potential, 

V.c{f) = -A-^—, (64) 
K r 

is of order A, and hence is an irrelevant perturbation in the RG sense that it vanishes as A — > 0. This is the choice made 
in the KSW scheme. In contrast, if we treat A^r as a low-energy scale and express it in units of A, writing Att = AAjr, 
then the rescaled potential is independent of A and so forms part of any fixed point of the RG. In Weinberg's power 
counting jj, this choice means that the potential (in momentum space) is of order It must thus be iterated 

to all orders when solving the Schrodinger equation, along with the leading contact interaction, and so this choice 
corresponds to the WvK scheme. 

In the ^Sq channel, the DW's of the Yukawa potential tend to constants as i? — > and so have the short-distance 
form 

^l{p,R) ^ I^{aT,/p,m^/p). (65) 



This corresponds to setting a — 1 and F{XR) = 1 in Eq. H59|l above. The resulting RG equation for Vs is given by 
Eq. 161|) . with a — I and three low-energy scales, p, a-^ and niT^. 

The expansion of the potential around the trivial fixed point has the form 

Vs^Y. CuranAPmfa'^p^^. (66) 

k.ni,n 

The RG eigenvalues of these terms are p — 2fc + m + 2n + 1 where fc, m and n are non-negative integers. The 
corresponding power counting assigns them orders Q'^ where d = p — 1 \2Ci l2l| . As in the case of a pure short-range 



° This is just the coordinate-space version of the rescaling discussed above. 

^ An alternative way to arrive at this result is to assign the nucleon mass an order l/Q in the power counting [II l4fill . This also leads to 
being identified as a low-energy scale of order Q. 
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potential, their eigenvalues start at p = 1 and so they are all irrelevant. This would be the appropriate power- counting 
if all short-range interactions in the ^Sq channel were weak. 

In the RG framework, the WvK treatment of this channel corresponds to an expansion of the short-distance potential 
around the nontrivial fixed point Vso- This can be written in the form 

- - E C,„„A''mf (67) 

Vso k,m,n 

The RG eigenvalues of these terms are p ~ 2k + m + 2n — 1 where fc, m and n are non-ne gati ve integers. This is 
similar to the expansion around the nontrivial fixed point for a pure short-range potential l2fl|| in that the leading 
perturbation is a relevant one, with eigenvalue p = —1. The terms are in one-to-one correspondence with the terms 
in a DW or "modified" effective-range expansion 0,E3jEEE3- In it, all rapid energy dependence associated with 
the low-energy scales of OPE is factored out, to leave an amplitude whose energy dependence is controlled only by 
scales from the short-range physics. 

In this channel we therefore have a choice between the two schemes, both of which lead to consistent expansions 
of the low-energy physics. The KSW scheme suffers from poor convergence because of the smallness of Ajr, while the 
WvK one, by treating A^r as a low-energy scale, converges better but lacks a clear connection with ChPT. The most 
immediate signal of the latter problem is the contact interaction needed to renormalise the logarithmic divergence 
produced by the 1/r singularity of the potential. This is proportional to a^r, which contains two powers of in the 
chiral expansion but which must be treated as a single small scale in the the WvK scheme |2l|. This is discussed 
by Beane et al. T^, who conclude that the KSW scheme should be used in the ^5*0 channel. However, it is worth 
stressing that the problem is a lack of consistency with the chiral power counting for other effective operators, not an 
internal inconsistency of the DW effective-range expansion. 



B. Tensor OPE 



The nonperturbative nature of the tensor piece of OPE at short distances means that the RG analysis is rather 
different from that for the central piece just discussed. In Sec. II we have seen that the wave functions in all partial 
waves (and for both attractive and repulsive potentials) have the same power-law dependence on x ~ pr at short 
distances. If we choose R to be small enough that it lies in the asymptotic WKB region, then the uncoupled DW's 
can all be written in the form 



Mp. R) - ^(K/p) {pR)-^'^F{KR), 



(68) 



where F^Xt^R) is the non-power-law part of the wave function. In the isospin-triplet waves, where the tensor potential 
is repulsive, this function is 



F{X„R) = exp 



-2 



X-^R 



(69) 



In the isospin-singlet channels, we have 



F{X„R) = cos 



X-^R 



+ 7 



(70) 



where 7 is an energy-independent phase. Here, for simplicity, I will continue to consider the chiral limit of the tensor 
OPE potential. Allowing for a finite will not alter the form of the short-distance wave functions, but it will 
introduce a dependence on rriT^/p into the normalisation constants A/". 

For the coupled channels, we have two independent solutions which satisfy the two boundary conditions as needed 
to generate an orthognal set of DW's, namely that as r ^ the exponential piece should be regular and the oscillatory 
piece should have an energy-independent phase 7. These solutions can be written 



^PL^{P,R)-^^+{K|p) {pR)-^'^F+{X^R)e+ + MUKIp) {pR)-^'^F_{X^R) b_ 
where the short-distance forms of the wave functions are 



1,2, 



(71) 



i^+(A^i?) =exp 



-2j^ 



F^{X^R) = cos 



2J — +7 

X 



(72) 
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and K± and e± are the eigenvalues and eigenvectors of the potential matrix K^ , defined in Eqs. (|54I55|I . 
Matrix elements of the short-range potential between DW's can be written 

{^L^\Vs\^Lk) ^ J2 ^f'UKIpTFo,{KR)*Vso.l3{p,K,KR)Ffs{KR)Nkp{K/p), (73) 

where the strength has been expressed as a matrix 'V s{p, A^r, A, R) using the basis of the eigenvectors e±. 
By analogy with the single-channel case, I define a rescaled potential V5 as a matrix with elements 

= ^ (Ai?)-i/2 F^{KR)* VsMp, a., a, R)F^{X^R). (74) 

It is also convenient to define the matrix N(At/p) with elements 

Nai=^f^a{K/p). (75) 

Using this in the two-component version of Eq. (|58|l leads to the RG equation for V5, 

^l?/=^^ + ^-l^+9^^ + N(A.)N(A.)t V^,. (76) 

d\ dp 9A^ 2 1 - 

Apart from the additional matrix structure, this equation has the same form as the RG equation (|61|l that governs 
the potential in the uncoupled channels. In particular, the coefficient cr, which arises from the power-law behaviour 
of the wave functions at short distances, is the same in all cases (attractive, repulsive and coupled). This means 
that the scaling behaviours of perturbations are the same in all cases, although obviously the forms of any nontrivial 
fixed points are different. In writing down these equations, I have ignored the possibility of bound states in the 
attractive channels. I discuss below the modifications needed to take account of these, but they do not affect the 
scaling behaviour. 

As always, each of these RG equations has a trivial fixed-point solution Vs — 0. The expansion around this has the 
form 

m,n 

where n and m are non- negative integers, and the RG eigenvalues of the terms are p = to + 2n 4- i. They are 
all irrelevant perturbations, in the sense that their eigenvalues are positive, and so the fixed point is stable. The 
appearance of noninteger anomalous dimensions should not be surprising in the context of the RG. The version 
developed in Ref. [23| shows that the scaling is controlled by the power-law behaviour of the wave functions at short 
distances and in general this is noninteger for potentials with an inverse-power-law form.^" In the present case, 
half-integer values appear because ~ (pi?)~^/^. 

For comparison, the terms in a pure short-range potential describing weak scattering in an S'-wave have RG 
eigenvalues p — 2n+ \ 20] . From this we see that the corresponding terms in the presence of a 1 /r'^ potential have 
eigenvalues that are smaller by subtraction of 1/2 and so they vanish more slowly as A ^ 0. This means that the 
effects of these interactions have been enhanced by the \/r^ potential. (In the language of the RG, they are more 
"relevant".) If we translate these results into the more usual power counting, we find that these terms have orders 
Qm+2n-i/2 ^ and so they have been promoted by half an order compared with the power counting one would get from 
naive dimensional analysis (Q™"'"'^"). 

In the S"-, P- and D-waves, the radii at which the OPE and centrifugal potentials become comparable are of the 
order of 1 fm. By choosing the factorisation radius R to lie significantly inside this, we can ensure that the short-range 
potential acts in a region where the wave functions have reached their asymptotic, WKB form for the \/r^ potential. 
In the '^Pi and '^£'2 channels, for example, this occurs for radii of the order of 0.1 fm or less. At these distances the 
wave function in the ^Pi channel is is highly suppressed by the repulsive \/r^ potential. Somewhat larger radii, of the 
order of 0.4 fm, can also be used since the wave functions still have energy-independent forms, at least for energies 



Indeed other examples of noninteger anomalous dimensions have recently been found in the context of three-body systems l5Ci . 
By neglecting the term on the LHS and changing variables to y = \J k/x, Eq. JHJ can be put into the form of Bessel's equation. This 
shows that the energy-independent solutions can be expressed in terms of ordcr-2/ + 1 Bessel functions of y. For large y, these tend to 
the asymptotic WKB expressions discussed above. 
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up to about 250 MeV, but even here the ^Pi wave function is already small, down by a factor of 5-10 compared to 
its size in the region 0.7-1 fm. 

The corresponding short-range potential must therefore be enhanced by a large numerical factor if we choose R 
less than 0.4 fm. However one should remember that the potential is not a physical quantity. In any observable, Vs 
always appears multiplied by two short-distance wave functions. Indeed the radius R is arbitrary and the form of 
the i?-dependence of Vs was chosen to ensure that scattering observables are independent of it. It is the size of the 
physical effects of a term in Vs that matters and, in particular, whether these require that it be iterated to all orders, 
or allow it to be treated perturbatively. Such questions are answered by the power counting that the RG analysis 
provides. 

Note that, despite also using a radial cutoff, the approach here is quite different from that Pavon Valderrama 
and Ruiz Arriola 0, 123L : those authors consider the limit as their radial cutoff tends to zero and so they set all 
irrelevant perturbations to zero. As a result, their predictions for scattering observables just depend on the long-range 
pion-exchange potential and a small number of short-distance parameters associated with any relevant or marginal 
terms. 

The results above apply to all partial waves where tensor OPE is treated nonperturbatively. Hence in waves with 
nonzero orbital angular momentum, where the leading terms are naively of order the orders of these terms are 
much lower than dimensional analysis would suggest. This agrees with the conclusion of Nogga, Timmermans and van 
Kolck |l4| that, based on their numerical analysis, short-range terms must be promoted in channels where the tensor 
potential is attractive. The RG analysis here makes quantitative the degree of promotion involved by determining 
the power counting for all terms in the double expansion in powers of energy (p^) and the coupling scale (Ajr). It also 
shows that the effect is present in repulsive as well as attractive channels. 

The scattering in the ^Si-^Di channel is strong at low energies and so the trivial fixed point is not an appropriate 
starting point. Instead we need to find a nontrivial fixed point and expand around it. This is most easily done by 
rewriting the RG equation as in Eq. (|62|l . The linear nature of this equation makes it straightforward to find the 
perturbations around the fixed point that scale with definite powers of A. The resulting expansion is 



where n and m are again non-negative integers. As in the case of a pure short-range potential |2flj, the nontrivial 
fixed point is unstable, with one negative eigenvalue. Terms in this e xpan sion can be related to the terms in a DW 
effective-range expansion, analogously to the examples studied in Refs. |2lLl23| . The RG eigenvalue for a general term 
in Eq. H78(l is p — m + 2n — 1/2, which should be compared with p = 2n — 1 for the pure short-range case [23|. This 
shows that the terms in the expansion here have been demoted by half an order (that is, they are less important) 
compared with the corresponding terms without the long-range potential. 

Although strictly the RG equation H61() only applies to channels where the tensor force is repulsive, the scaling 
behaviour is in fact the same for the attractive and coupled channels. The only difference is that that attractive 1/r^ 
potentials give rise to deeply bound states that lie outside the domain of our EFT. We should therefore cut them off 
ai E = — A^/M„, as in the case of the attractive inverse square potential This adds (5- function terms to the RG 
equation at the values of A where bound states fall outside the cut-off. These lead to step discontinuities in Vs at 
these points, which can be thought of as jumps to different branches of the fixed-point potential Vso- 

The existence of multiple branches of the potential is a consequence of the oscillatory nature of the short-distance 
wave functions. In order to make these well defined, we had to choose a particular self-adjoint extension of the long- 
distance Hamiltonian by fixing the phase of these oscillations. As in the 1/r^ case, the energy-independent short-range 
interaction has the effect of changing the self-adjoint extension |23|. However, in that example, the scale-free nature 
of the potential and the associated Efimov effect |^ (an infinite tower of geometrically spaced bound states) mean 
that the different choices lie on a limit cycle of the RG. As a result the leading short-range term forms a marginal 
perturbation which changes the starting point on that cycle. In contrast, only a discrete set of extensions of the 1/r^ 
Hamiltonian lead to scale-free systems with bound states at zero energy. The different branches of V50 correspond to 
this set of extensions. 

The use of the DW basis in this RG analysis makes it straightforward to expand the potential in terms of perturba- 
tions that scale with definite powers of A. In Ref. 14] the cut-off was applied to a plane-wave basis. Such a cut-off has 
two effects: it both regulates the short-distance interaction, and it removes the singularity of the long-range potential 
at the origin. The second aspect means that changing the cut-off has the effect of changing the self-adjoint extension 
that determines the long-range behaviour of the DW's. After fitting to low-energy scattering observables, the resulting 
short-distance potential displays dramatic oscillations due to the changing number of bound states of the long-range 
potential. Similar oscillations are also seen with the radial cut-off of Ref. [l^. However, with care, it should still be 
possible to determine the power counting using such cut-offs, as illustrated by the analysis of attractive three-body 




(78) 



systems in Ref. |53 |. 
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In spin-triplet channels without low-energy bound or virtual states (in other words, all except "^Si-^Di), we can 
define the DW basis by picking an initial extension that does not produce a low-energy bound or virtual state and then 
we can expand the short-range potential around the trivial fixed point. The terms in this expansion are all irrelevant. 
This implies that, provided we pick an initial extension that gives weak low-energy scattering, any dependence of 
the low-energy phase shifts on this choice must be small. Otherwise the leading irrelevant perturbation could not be 
equivalent to a change in the extension. 

This weak dependence of scattering observables on the choice of extension, at least well away from the ones that 
generate low-energy bound states, can be seen in the numerical results of Ref. ,14] . particularly in Figs. 10 and 12. 
These contain long "plateau regions" where the short-distance interactions are small and only weakly dependent on 
A. These values of the cut-off correspond to extensions that generate only deeply bound states, outside the scope of 
the EFT. In contrast there are also narrow ranges of A where the regulated tensor potential produces a low-energy 
bound or virtual state. To describe the observed weak scattering, the potential in these regions must be supplemented 
by the nontrivial fixed point and its relevant (unstable) perturbation, hence the very large counterterms needed. 

The results so far can be applied to the S-, P- and D- waves where, as we saw in Sec. II, the tensor OPE needs to be 
treated nonperturbatively. As discussed above, the wave functions in these channels attain their energy-independent 
short-distance forms for radii of 0.4 fm or less. These forms are controlled by the Ijr'^ tensor potential which is 
stronger than the centrifugal barrier for radii smaller than about 1 fm. In _F- waves and above, in contrast, the OPE 
and centrifugal potentials become comparable at radii of the order of 0.1-0.2 fm or less. The asymptotic forms of the 
waves controlled by the tensor potential are thus reached only for radii much less than 0.1 fm, far beyond the domain 
of validity of our EFT. Moreover, the amplitudes of the wave functions for such small radii will be strongly suppressed 
by the centrifugal barrier, at least for momenta of the order of m^.^^ In effect the strong centrifugal barrier in high 
partial waves "protects" low-energy waves from probing the nonperturbative region. 

We can take advantage of this by choosing our factorisation scale R for the high partial waves to be of the order 
of 0.4 fm, so that it lies in the region where the centrifugal potential dominates over both the tensor OPE and the 
on-shell energy. This allows us to treat OPE as a perturbation and to expand the DW's in Eq. 158|l in powers of 
In this region the waves still have the r' form produced by the centrifugal barrier and so scaling behaviour of the 
short-range potential is the same as that in the presence of the 1 /r^ centrifugal potential alone . In particular the 
resulting power counting is just the usual one given by naive dimensional analysis. Of course, at high enough energies, 
the waves will penetrate the barrier and this perturbative treatment will break down. However, in the high partial 
waves, the results in Sec. IIC show that this only happens at scales that lie outside the domain of our low-energy 
effective theory. 

IV. CONCLUSIONS 

The chiral limit of the tensor OPE potential has a 1/r^ form. In this work, I have studied it using techniques 
developed in atomic physics for solving the Schrodinger equation with inverse-power-law potentials piill7lll8l |. These 
lead to analytic solutions which are constructed as expansions in Bessel functions whose orders satisfy an eigenvalue 
equation. In each channel, there is a critical value of the product of the momentum and coupling strength above 
which these eigenvalues become complex. This sets a limit on the range of energies for which the tensor potential can 
be treated perturbatively. These values are the same for both attractive and repulsive l/r^ potentials. 

I have determined the critical dimensionless couplings for the tensor interaction in low- lying partial waves. In the 
^Si-^Di and ^Pq channels, the corresponding breakdown scales are ~ ttItt or less. In the other P- and _D-wavcs, the 
scales are of the order of 400 MeV. These results imply that, for the energies relevant to nuclear physics, OPE should 
be treated nonperturbatively in these channels. They explain why Fleming, Mehen and Stewart 12] found that the 
perturbative KSW treatment breaks down in these cases. In contrast, the scales for the higher partial waves all lie 
well above 1 GeV and so perturbation theory should be valid for them. 

In the context of the RG, the nonperturbative treatment of OPE can be justified if we identify the scale Att 
controlling its strength as a low-energy scale. The resulting RG analysis of the ^5*0 channel, where only the central 
Yukawa piece contributes, leads to power countings that are similar to those found for pure short-range interactions. 
There is a nontrivial fixed point which describes systems with strong scattering at low energies. The terms in the 
expansion around this have RG eigenvalues /9 = —l, 0,-1-1, ••• and so are of order Q"^ where d = /9— 1 = — 2,— 1,0,--- 



The only exception would be if we were to choose an extension of the attractive potential that leads to a low-energy bound or 

virtual state which is trapped inside the barrier. However there is no reason to do so in the context of nucleon-nucleon scattering in 
high partial waves. 



18 



[2l|. This is similar to the power counting for the expansion around the nontrivial fixed point for pure short-range 
forces 01 120( . In both cases the terms in the expansion correspond directly to terms in an effective-range expansion. 

Here I have used this RG method to study the scaling behaviour in the spin-triplet channels. This is controlled by 
the power-law dependence of the DW's near the origin which, in turn, follows from the singularity of the long-range 
potential at the origin. It is the same for attractive and repulsive potentials. In the S-, P- and D-wave channels, 
where the tensor OPE force needs to be treated nonperturbatively, the expansion around the trivial fixed point leads 
to RG eigenvalues p = ^jfjfr''j corresponding to orders Q"^ with (i = — These are promoted by half an 
order compared with naive dimensional analysis in an S'-wave and by many more orders in higher partial waves. This 
provides a quantitative measure of the effect observed in Ref. 0| . 

In the channel we need to expand around the nontrivial fixed point. The corresponding RG eigenvalues 

are p=— - and so there is one relevant perturbation. This is similar to the pure short-range case, except 
that the terms in the expansion have been demoted by half an order. 

Finally, in the higher spin-triplet waves, one can treat the tensor OPE potential perturbatively at low energies. 
The corresponding scaling behaviour is determined by the region where the centrifugal barrier dominates, and so the 
power counting is just that given by naive dimensional analysis. 

These results show that it is possible to set up a consistent EFT embodying the WvK scheme where OPE is treated 
nonperturbatively in low partial waves. However, as already remarked, a central element of this is the identification 
of Att as a low-energy scale. Since this scale is built out of quantities that are treated as high-energy scales in ChPT, 
this analysis leaves open the question of how to make this theory consistent with chiral expansions of other effective 
operators, such as those for EM or weak couplings. Further work is needed to address this. 
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